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ABSTRACT 

We model solar coronal mass ejections (CMEs) as expanding force-free mag- 
netic structures and find the self-similar dynamics of configurations with spatially 
constant a, where J = «B, in spherical and cylindrical geometries, expanding 
spheromaks and Lundquist fields respectively. The field structures remain force- 
free, under the conventional non-relativistic assumption that the dynamical ef- 
fects of the inductive electric fields can be neglected. While keeping the internal 
magnetic field structure of the stationary solutions, expansion leads to compli- 
cated internal velocities and rotation, caused by inductive electric fields. The 
structure depends only on overall radius R{t) and rate of expansion R{t) mea- 
sured at a given moment, and thus is applicable to arbitrary expansion laws. 
In case of cylindrical Lundquist fields, magnetic flux conservation requires that 
both axial and radial expansion proceed with equal rates. In accordance with 
observations, the model predicts that the maximum magnetic field is reached 
before the spacecraft reaches the geometric center of a CME. 



Introduction 



Expansion of magnetic clouds is one of the basic problems in space physics and geo- 
physics, related to propagation of solar disturbances, the coronal mass ejections (CMEs) 



through interplanetary medium (Forbes et al. 2006). As the space craft passes through the 



cloud, the magnetic field strength is higher than average, the density is lower, the magnetic 
pressure inside the cloud greatly exceeds the ion thermal pressure and the magnetic field 



direction changes through the cloud {e.g. Burlaga et al. 1982). 



There are two basic models of magnetic clouds: (i) magnetic flux ropes, still connected 



to the Sun at the distances of the Earth (Burlaga et al. 1981 Cargill et al. 1995 Chen and 



Garren|1993 ); (ii) disconnected entities of spherical topology, which in the simplest paradigm 
can be modeled as spheromaks {e.g. 



Ivanov and Kharshiladze 1985: Vandas, Fischer, and 



Geranios|1991 ). These are stationary solutions, while the CMEs expand while propagating in 
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the Solar wind. Previous attempts to find a structure of expanding spherical configurations 
have led to breakdown of force-free condition through appearance of pressure and inertial 
effects QLow||1982t |Ivanov et qfll989| |Vandas et aE]|1997[ ). In this paper we find the structure 
of expanding magnetic configurations that remain force-free. In the force-free approach the 
dynamics of the system is dominated by the magnetic field and any plasma pressure or 
inertia terms are negligible, for that reason we shall not include in this study the equations 
of continuity of the fluid and we shall focus only on the magnetic field and on the induced 
electric field. 



The importance of spheromaks in spherical geometry ( Chandrasekhar and Kendall 1957) 



and of Lundquist (1951) fields (in cylindrical geometry) as examples of clouds with self- 



contained magnetic field stems from the fact that these solutions are, in some sense, the 



minimal energy states and thus are stable to many ideal and resistive instabilities. Woltjer 



(1958) and Taylor (1974) formulated the plasma relaxation principle, according to which 



plasma evolution conserves magnetic helicity and, if the magnetic field is the only important 
quantity, it evolves to a force-free state with constant a, where J = «B, where a has some 



similarities to thermodynamic temperature, so a relaxed state should have constant a (Bellan 
20061). 



2. Expanding Spheromaks 

2.1. Magnetic Field of a Spheromak 

Spheromaks are stationary force-free configurations of plasma satisfying the condition 
J = aB with spatially constant a. They are solutions of Grad-Shafranov equation in spher- 
ical coordinates when the poloidal current is a linear function of the magnetic flux. Below 
we will use "the basic spheromak" solution, corresponding to axially symmetric dipolar-like 
fields (Chandrasekhar and Kendall 1957). 

Given the magnetic flux function P(r,9) and poloidal current function I(P), the mag- 
netic field is 

VP 21 (P) 

B = ^ xe ^ + ^ 1 
r sin 6 r sin 6 

When only magnetic force is present, the magnetic flux function satisfies the Grad-Shafranov 
equation 

A*P + All' = 0, (2) 

where 

a- = a; + ^* A= (3) 

r z sin v 
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is the Grad-Shafranov operator. For linear dependence / = aP/2 with spatially constant a, 
Equation ^ has separable solutions in term of spherical Bessel functions j n (ar). The first 
order solutions are then 



P = B -ji(ar) sin 2 9, 
a 

B r = 2Bq— cos 8, 
ar 

B e = -B Jl + arj[ sine, 
ar 

B^ = B o3l sin 9. (4) 

The first order spherical Bessel function j\ can be expressed in terms of elementary functions, 
ji(x) = sinx/x 2 — cosx/x. Parameter a is related to the size of a spheromak, defined by the 
surface where radial magnetic field is zero, solution of j\ = 0, R = C a /a, C a = 4.49. 



2.2. Radial Expansion of a Spheromak 

We now study the self-similar, non-relativistic expansion of axially symmetric force- 
free structures with spatially constant a-parameter. Below we discuss a procedure that 
allows to derive non-relativistic, time- dependent force-free solutions from a given stationary 
force-free configuration. First, we note that parameter a has dimension of inverse length, 
and thus is related to overall size of the magnetic structure R. For expanding solutions 
we take a = a(t), but still spatially independent. According to self-similar prescription, 
we take all radial dependence to be in a form oc ra oc r/R(t). The induction equation, 
V x E + d t B = 0, as well as magnetic magnetic flux conservation, require that in a time- 
dependent case P = (a/«o) 2 -Po) where «o is the value at some initial time and Pq is a solution 
of a stationary force-free problem (e.g. given by Equation (|4])). Next, we are looking for 
ideal solutions, where electric field is always perpendicular to magnetic field and thus can be 
written as E = —V x B. Since we are looking for radially self-similar solutions, V = Ve r (cf. 



Low 1982, and Equation (|7j)). This is the velocity of radial expansion. The drift velocity of 
the magnetic field lines is given by v = E x ~Q/B 2 and is by definition perpendicular to the 
magnetic field lines. These two approaches give different values for V and v. Using vector 
identities we can conclude that the difference between two velocities is a vector parallel to 
the magnetic field lines. Indeed 

ExBVxB VB 
-= — = -— xB = V-^B. (5) 

There is a freedom of choice between these velocities without altering at all the force-free 
solution. Equilibria of this type, which contain also some matter however, do not allow 
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transition from one velocity to the other without changing the dynamics of the problems just 
by adding a velocity component parallel to the magnetic field. In that case, an azimuthal 
component in the velocity will give rise to inertial terms in the force balance equation. In 
this study we focus in force-free equilibria where the dynamics of the system are dominated 
by the magnetic field and inertia is negligible. Choosing self-similar scaling, V = —rdt In a, 
we find for axially symmetric magnetic field, given by Equation 0, the following electric 
field 

(X 

E = r-e r x B, (6) 
a 

where dot denotes differentiation with time. It follows from Equation ^ that: 

B(B-e r )-e r .i? 2 
v = — rd t \na. (7) 

In deriving the solution we neglected dynamical effects of electric field, in particular dis- 
placement current and electric charge density, so that the force-free condition involves only 
magnetic field and electric current, J x B = 0. This is indeed a reasonable assumption for 
the solar system where the velocities are non-relativistic. 

For the basic spheromak solutions Q, the electromagnetic fields become 

B r = 2B a-^- cos 6, 

R r h + ar fi ■ a 
Bg = —B a 5 sinS, 

B6 = B j 1 ( — ] sin 9, 



a J 



E r = 0, 



_, __. aa a 
E e = -So^-rji sinf = -B^-r, 
chq a 

P r> ■ h + arj[ a 

Es = —B a k sm& = B e —r, 

Q!q a 



E = ( -r ) e,. x B 



We stress that the fields remain force-free, J x B = 0. 



2.3. Velocity Structure of the Solution 



Expanding force-free spheromaks have complicated internal velocity structure. At each 
point there is a well defined electromagnetic velocity v = E x B/Z? 2 , normalised to c, 
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Equation ([7]). The solution can be parametrized by the velocity of the expansion of the 



boundary, which is v g = = and v r = R(t)d t ln(l/a) = v , see Figures lpp 



o 

> 

> 




r/R 



Fig. 1. — Radial velocity v r , Vg, as functions of r for different 9. In the equatorial plane 
radial velocity increases linearly with radius, v r (9 = it/ 2) = (r/R)vo, while close to the axis 
velocity remains small in the bulk, sharply increasing to vq near the surface. 



Since ^ 0, the expansion of a spheromak induces rotation, with the different hemi- 
spheres rotating in the opposite direction (Figure [4]). The maximum toroidal velocity is 
0.33t>o reached at 9 = 0.53 and r/R = 0.71. For r <^ R the angular velocity of rotation is 
co ~ ±C a v o sin 29 The flow lines are plotted in Figure [5| Close to the symmetry axis, 
velocity at small radii is directed away from the axis.. Total velocity at small r R is 
v ~ Vqt sin9/R. The flow lines are orthogonal to the surface, matching the overall expansion 
of the spheromak. 

Using expressions (J8l) for electromagnetic fields, we can calculate the total energy of the 
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0.0 0.2 0.4 0.6 0.8 1.0 

r/R 

Fig. 2. — Velocity v$ as a functions of r for different 9. Flows are antisymmetric with respect 
to the equatorial plane, vq becomes zero at r = 0.61-R (root of the equation E$ = 0) 

expanding spheromak 

E tot = j dV 

where C\ = 2.5 x 1CT 3 and C2 = 8.7 x 10~ 4 are combinations of power laws and trigonometric 
functions of C a . The second term proportional to velocity squared is the contribution of the 
electric field and can be ignored for non-relativistic expansion. 

It may be verified that the total magnetic helicity, H = J dV~B ■ A = 1.4 x 10~ 2 BqRq, 
as well as the toroidal magnetic flux J 7 = j B^rdr = 5.26BqRI are independent of time. At 
the same time, the total electromagnetic energy decreases during expansion oc 1/R. This 
highlights an important point, that it is impossible to conserve both electromagnetic energy 
and magnetic helicity. In case of ideal expansion the energy decrease is due to pdV work 
done on the surrounding plasma. Magnetic pressure on the surface cx 1/-R 4 , energy oc 1/R. 

We can compare the self-similar spheromak solutions to the internal energy of adiabati- 
cally expanding polytropic gas with p oc p 7 , E oc pV oc R?^~^. Thus, effectively, spheromak 
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r/R 

Fig. 3. — Velocities as functions of r for different Flows are antisymmetric with respect 
to the equatorial plane. 




0.0 0.2 0.4 0.6 0.8 1.0 ' 0.0 0.5 1.0 1.5 

r/R 6 

Fig. 4. — Angular velocity of rotation Vt = v^/{r sin 0) as function of radius and polar angle. 



behaves as a polytropic gas with 7 = 4/3. This is related to the results of Low (1982), who 
found self-similar MHD solutions for expanding cloud if the fluid obeys polytropic law with 
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Fig. 5. — Plasma flow lines. On the surface flow lines become radial to match the overall 
radial expansion. The stationary point is 9 = 0, r = 0.61-R (root of the equation = 0). 



7 = 4/3. In this case both magnetic and thermal components of plasma behave similarly 
( Gourgouliatos and Vlahakis 2010). 
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3. Toroidal Configurations 



In the case where magnetic clouds are connected to the Sun (Chen and Garren 1993), 



their topology shall resemble more an arcade, and their mathematical description is more 
accurately approximated by a force-free magnetic field confined in a torus. A number of 



studies (Burlaga et al. 1981 Burlaga 1988; Lepping, Burlaga, and Jones 1990) have con 



sidered solutions in cylindrical geometry as an approximation to the toroidal structures of 
magnetic cloud. 

In this section we examine the effect of expansion in this type of solution by applying 
self-similarity. Let us assume a torus (Figure [6]) of large radius Rq and small ro, the centre 
corresponding to the large radius is O and the centre corresponding to the small radius lies on 
the circle (O, Rq) and shall be called O'. Let us attach an orthogonal system of coordinates 
to 0'\ this coordinate system resembles a cylindrical coordinate system with the difference 
that the role of the axis z is taken by the circle (O, Rq), the coordinates in this system are 
(zu,(j),z), where < w < ro, < < 2tt and —nRo < z < ttRq; the field inside the torus 
has components B mi and B z . Such fields have been studied by Miller and Turner (1981), 



who considered force-free structures in toroidal coordinates and found that for large aspect 



ratios Rn/r 3> 1 the solution reduces to Lundquist (1951) fields. The first few terms of the 



Ro 



Fig. 6. — The geometry of the torus, the large radius Rq writes a circle with centre O and the 
small one ro a circle with centre O' . We attach the orthogonal coordinate system (w, <p, z) 
to the circle (O r , r ). 

Miller- Turner fields are 

B^.t = -=-5- ( - J (aw) + — ) sin 9 , 
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Ji(aw) \- (j (aw) - r ^) cos6> 

alio V aw J 

B z ,t = B Q [J (aw) - -^-Fcos6» 



Brh.T = Bn 



o 

F=—J (aw) + - — Maw). 10 

2r 2Ji(ar ) 

Let us now consider a thin section of the torus which can be approximated with a cylinder. 
In the large aspect ratios limit, Ro/r ^> 1 we only keep the Lundquist field. 

B w ,s = , 



B<f>,s = B J i (aw 

B ZjS = B J (aw), (11) 



where Jo,i is the Bessel function of the zeroth and first order. 

Magnetic clouds expand as they propagate away from the Sun. The simultaneous ex- 
pansion and propagation lead to an increase of both radii of the torus. The expansions of the 
small and large radii need to scale with the same factor to conserve the magnetic flux that 
crosses the equator of the torus and the magnetic flux that crosses a section of the torus: 
$i cx B^tqRq and <3> 2 oc B z ^Tq, so if R is multiplied with some factor so does r . 

This field shall experience the effect of expansion both in w and z directions. Let us 
consider an expansion velocity V = —w-e^ — z-e z . 



Similar to Section [272| the magnetic flux conservation of the expanding structure requires 

that 



/ a N 2 
B^e — B \ — 

/ a n 2 



B m ^ = 
Ji (aw \ 



B Z ^ = B [ — ) J [aw). (12) 



Then the electric field is E = —V x B, which gives 



aa 

E w = —B —^-J 1 [aw)z. 
at 



Es = -B ^ J (aw) w , 
aa ( \ 

E z = B — =- Ji I aw ) w . (13) 
ai V / 
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The whole structure is force-free satisfying (V x B) x B = 0, indeed 



V x B = Bo^J^aw)^^ + B ^J (aw)e z = aB . 



(14) 



In addition, the divergence of the magnetic field is zero V • B = and the magnetic flux is 
conserved as it satisfies Faraday's law of induction V x E + <9 t B=0, indeed 

OLOL 

V x E = — B — [(Ji(azu) + awJ (aw))e ( f > — (awJi(aw) — 2J (azu))e z ] = — <9 t B . (15) 



We stress that though we focus in on a thin slice of the torus area near z — 0, where 
the V z component vanishes, we do take into account the fact that both the small and the 
large radii of the torus are increasing. The latter gives an expansion in the z direction which 
creates a minor electric field in the w direction which is necessary to conserve the magnetic 



Again we evaluate drift velocity of the field lines v = E x B/J? 2 , whose components are 



The components of v shall match V at the boundary. Indeed, the v w component is mathced, 
then from the (ft component we find that the cylinder shall end either at a root of J\ or J 
and then from the z component we find the acceptable solution is that of a root of J , thus 
the end of the rope is defined by the position where aw equals a root of Jo. 

Formally, the (ft and z velocity components are proportional to z and this breaks the 
translational symmetry of the Lundquist field. This is the outcome of choosing an arbitrary 
cross section and setting the velocity profile proportional to the distance from there. The 
physical system actually consists of a torus and as a result of the expansion an observer on a 
given point in the torus will see the other points to move away; and their distance measured 
on the z coordinate will increase proportional to z. Given the arbitrary choice of the point 
z = we have the freedom to set this point at the place where the detector crosses the 
magnetic cloud, provided that its trajectory is perpendicular to the axis. The detector will 



flux. 



a 

v w = w , 

a 




(16) 
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see a structure where the values are given for z — 0, the magnetic field is that of Equation 



(12), the electric field though is 



Ed>,o = —-do — r Jo a^ 7 w > 
"o V ; 
cm / \ 

■S*,o = Bo— * J\ (aw U7 . (17) 
"o V ; 



Finally the velocity of the fluid is 



a 

Vzufl = 



= ^ 2 ,o = . (18) 
We remark that similar results for the cylindrical Lundquist type fields have been obtained 



recently and independently by Dalakishvili et al. (2010). 



4. Matching to Driven Shock Solutions 



The magnetic clouds expand while propagating away from the Sun (Klein and Burlaga 



1982). Next we consider how the expansion rate scales with the distance from the Sun. 
This depends both on the structure of the Solar wind and on the conservation or not of the 
magnetic flux and helicity of the cloud, during its evolution, which in turn is related to the 
connection or not of the cloud to the Sun. A spheromak-type cloud is a disconnected entity 
which ideally conserves its magnetic flux and helicity, however as it expands its electromag- 
netic energy cannot remain constant. Magnetic flux tubes, however, are attached to the Sun, 
so injection of magnetic flux or helicity is more likely. 

Let the spheromak have magnetic flux $ enclosed within a volume of size R; the magnetic 
flux has dimensions M 1//2 L 3//2 T _1 . The magnetic field on the surface is not constant, but for 
simplicity we assume that it is of order Bq = $/i? 2 . Let the external medium have density 
p C xt, which has dimensions ML~ 3 ; the other variable is time whose dimension obviously is 
T. From dimensional analysis we find that expansion of a spheromak obeys a law: 

R - Pcxff 1 ^ 1 / 3 • (19) 

A somewhat different approach leading to the same result comes from balancing ram pressure 
of a strong forward shock, p ext (dR/dt) 2 , with magnetic field energy density on the surface 
B 2 ~ ($/i? 2 ) 2 . 
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For constant mangetic flux $, and assuming that external density scales with radius as 
p cxt oc R~ a , the self-similar expansion requires 

R oc t 2/(6 ~ a) (20) 

In a constant density medium this gives R oc t 1 ^ 3 . Note, that for self-similar scaling the 
energy in shocked plasma, oc R 2 p ext R 2 oc £ _2 /( 6_a ) decreases with time, as t -1 / 3 in constant 



density, (cf. Farrugia, Osherovich, and Burlaga 1995, Equation (44)). This can be achieved 



e.g., if the shocked plasma cools radiatively, but the energy loss should be tuned to the 
self-similar dynamics, which is not justified. 

Thus, generally, if a magnetic bubble drives a non-radiative shock, the surface pressure 
inside the bubble falls off faster with time than the post-shock pressure and the self-similar 
assumption breaks down. At every moment the size of the bubble would adjust to the 
post-shock pressure in a non-self-similar way. For example, if initially the bubble is over 
pressurized with respect to a surrounding medium of constant density, expansion will drive 
a shock with r s oc t 2 ^ 5 ( Sedov||l959 ). The post-shock pressure scales as p oc t _6//5 and the size 



of the magnetic bubble R oc t 3//10 . Thus, the whole system of magnetic bubble plus external 
shock medium evolves in a non-self-similar manner, with the relative size of the magnetic 
bubble decreasing R/r s oc t -1 / 10 . Most energy is transferred to the outside plasma at initial 
stages of expansion. 



4.1. Expansion in a Wind Environment 



Motion of a CME in a wind is composed of translational motion and overall expansion. 
The translational velocities of the clouds vary in the range of 300km s _1 to 1000km s _1 with 
average velocity 478km s -1 ( Bothmer and Schwenn|[l998 ). This velocity is of the same order 
of magnitude as the Solar wind velocity, [j The local density at the cloud's location is then 
Pext = Po r o/ r CME5 where tcme ~ ^sw^- Substituting in Equation (19) we find that 



R 
R 



Po 



Po 



-1/6 -1/3 



r, 



o u sw - 

UQ i/ 3 $l/3 r 2/3_ 



(21) 



Bothmer and Schwenn (1998) have found R oc r 



„0.78±0.10 



close to the value of 2/3 predicted 



by Equation (21 ) 



1 We remark that if we take into account the buoyancy force and the drag due to the motion relative the 
solar wind we find that the limiting velocity is of the order of the sound speed of the solar wind, which is a 
few 10 2 km s _1 depending on the temperature which is consistent with the numbers mentioned above. 
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5. Observational Signature 
5.1. Magnetic Field Measured 



The magnetic field measured by a magnetometer as a function of time depends on the 
relative velocity of the detector and spheromak boundary. Generally, the velocity of the 
detector with respect to the spheromak's center is dominated by the advection velocity of 
the spheromak with the Solar wind, which is typically a factor of several larger than velocity 
of expansion of the spheromak in the winds frame. 

Observations of CMEs often indicate that the maximum of B-field is reached before the 



spacecraft reaches the geometrical center (Burlaga 1995 Farrugia et al. 1993). It was long 



suspected that is a consequence of the expansion of the magnetic cloud while it moves past 



the spacecraft (Burlaga 1995, Section 6.5.2): when the detector reaches the middle of the 



spheromak, the overall normalization of the magnetic field decrease due to expansion. Our 
model provides a quantitative description of this effect, see Figures [7] and [8} 



5.2. Plasma Velocity 

Observations suggest that the general trend of the velocity measured by a detector 



crossing the magnetic cloud, is to start from a large value and then decrease (Bothmer and 



Schwenn 1998) to a smaller but comparable velocity. This is consistent with the image of 



a magnetic cloud that propagates from the Sun and at the same time expands. We have 
assumed that the propagation velocity is constant while the expansion is self-similar R oc t n , 
see Figure |9j In general, clouds differ in the details of the velocity profiles, but most of them 
give measurements that are in agreement with this model. 



Discussion 



In this paper we discuss the self-similar solutions for expanding force-free magnetic 
structures in spherical and cylindrical geometries. Under the assumption of non-relativistic 
MHD, when the dynamical effects of the resulting induced electric field can be neglected, 
the structure remains force- free during expansion, J = qB. In case of relativistic expansion, 



when electric field has important dynamical effects (Prendergast 2005; Gourgouliatos and 



Lynden-Bell 2008), spheromak- type self-similar solutions may be found for constant expan- 



sion velocity (R. D. Blandford, private communication). In contrast, our solutions depend 
on R and R measured at a given time, and thus are applicable to any R(t). 
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Fig. 7. — Left Panel: Sections of a spheromak with a meridional plane. The trajectories 
of the detector crossing the spheromak are the dashed lines. Right Panel: Magnetic field 
measured by a detector flying through an expanding spheromak. The detector enters the 
cloud at t — 0, at t — 1 it crosses the axis: curves (1), (2); or the equator: curves (3), 
(4). We investigate four cases corresponding to the numbers of the trajectories shown in left 
panel. The solid curve (1) corresponds to a detector flying on the equator of the spheromak 
and passing from the centre at t — 1. The dashed curve (2) corresponds to a detector that 
flies parallel to the equator and crosses the axis at t — 1. The dashed-dotted (3) curve 
corresponds to a detector that flies parallel to the axis and crosses the equator at t — 1. 
Finally the dotted curve (4) corresponds to a detector flying along the axis and passing from 
the centre at t — 1. In all cases the maximum magnetic field is measured at t < 1 before the 
detector crosses the middle of the cloud. 



In case of spheromaks, the internal velocity develops azimuthal component, ^ 0, 
even though the surface expands purely radially. The fact that the magnetic field structure 



remains force-free during expansion is qualitatively different from the solutions of Low ( 1982 ) ; 



Farrugia, Osherovich, and Burlaga (1995), which were limited to pure radial motion. As a 



result, even if at the initial moment magnetic field is force-free, during expansion non-zero 
pressure forces appeared. 

We expect that the solution is dynamically and resistively stable, since stationary force- 



free solutions with constant a are very stable configurations (Bellan 2006), both to ideal 



and resistive instabilities, like tearing mode. Stability of expanding force-free structures 
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12 3 

t 

Fig. 8. — Magnetic field measured by a detector flying through an expanding cylinder. The 
solid curve (1) corresponds to a detector that flies and passes from the centre, whereas the 
dashed curve (2) correspond to a detector that enters the structure at an impact angle of 
7r/6 and misses the centre by half radius. 

should be investigated independently, but we expect them to be stable to global modes, since 
expansion-induced velocity shear will contribute to global stability (but may still result in 
development of local shearing instabilities). 

Expansion of magnetic clouds is qualitative different from the self-similar fluid outflows, 
since in this case we have an additional constraint, the conservation of magnetic helicity. 
Since the conservation of magnetic helicity and magnetic energy are, generally, not com- 
patible with each other, the self-similar expansion of spheromak can be connected only to 
one particular self-similar (Sedov-type) solution, with r oc t 1 / 3 . This solution has energy 
decreasing with time and thus requires strong energy losses. 

The model predicts that that the maximum magnetic field is reached before the detector 
reaches geometrical center of a CME, while the plasma velocity measured by the spacecraft 
smoothly decreases with time. Both these facts are in general agreement with observations. 
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Fig. 9. — Upper panel: The velocity measured by a detector that encounters an expanding 
magnetic cloud and crosses its axis. These results are valid for both a spheromak cloud 
provided the detector's trajectory is on the equator and for a cylindrical one provided the 
trajectory is normal to the axis. We have assumed that the cloud propagates at constant 
velocity and expands self-similarly so that R oc t n . The detector measures the combined 
effect. The detector enters the cloud at t = and exits at t = 1 where the velocity is 0.7 
of the initial. The solid line corresponds to a self-similar expansion with index n = 2/3, 
the dashed to n = 1/2 and the dotted to n = 1/3. Lower panel: The velocity measured by 



the Helios of the event on days 362-365 1978 (Bothmer and Schwenn 1998). For n < 1 the 



velocity profile in the cloud is consistent with the profile of the expanding model. 
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